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Abstract 

The effective action for local composite operators in QED3 is considered. The 
effective potential is calculated in leading order in 1/Nt (Nf is the number 
of fermion flavors) and used to describe the features of the phase transition 
at Nf = N CT , 3 < N CT < 5. It is shown that this continuous phase transition 
satisfies the criteria of the conformal phase transition, considered recently in 
the literature. In particular, there is an abrupt change of the spectrum of light 
excitations at the critical point, although the phase transition is continuous, 
and the structure of the equation for the divergence of the dilatation cur- 
rent is essentially different in the symmetric and nonsymmetric phases. The 
connection of this dynamics with the dynamics in QCD<± is briefly discussed. 
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I. INTRODUCTION 



Noncompact quantum electrodynamics in 2 + 1 dimensions (QED 3 ) with Nf flavors of 
four- component fermions |L],|2| is a gauge theory with rich dynamics, reminiscent of four- 
dimensional quantum chromodynamics (QCD4). 

Studying the Schwinger- Dyson (SD) equation for the fermion self-energy in leading order 
in 1/Nf expansion || showed the existence of a critical number of fermion flavors N cr , 
3 < N CT < 5, below which there is the dynamical breakdown of the flavor U(2Nf) symmetry 
and fermions acquire a dynamical mass. These conclusions were confirmed in the lattice 
computer simulations of noncompact lattice QED 3 [Q. 

The presence of a critical N CI in QED 3 is intriguing especially because of a possibility 
of the existence of an analogous critical Nf = N cr in QCD4, suggested by both analytical 
studies |5|-|8| and lattice computer simulations of QCD^ |p|-|lT|. 

However there is still controversy concerning the existence of a finite N CT in QED 3 : 
some authors argue that the generation of a fermion mass occurs at all values of Nf in 



QED 3 [|l2| , [l3l . Despite quite intensive studies of this issue |Tj,|T5[, the problem is still 



open. The main difficulty is that at present there is no systematic approach to studying 
nonperturbative strong coupling dynamics in gauge theories. A possibility to get insight 
into the phase transition at Nf = N CT might be a low energy effective action approach, 
which has been very successful in studying the phase transition at finite temperatures in 
QCD^ [|TBJ. Such an approach is useful when there are few relevant (massless) degrees of 



freedom at the critical point. For example, at T = T cr in QCD4, the relevant degrees of 
freedom are Nj (or Nj — 1) pseudoscalar and Nj scalar mesons: because of a non-zero 
temperature, and therefore non-zero Matsubara frequences for quarks, all other excitations 
are effectively massive at T = T cr . As a result, the effective theory is given by a three 
dimensional renormalizable lagrangian with the complex matrix field (pi, i, j = 1, 2, . . . , Nf, 
describing 2Nj scalar and pseudoscalar mesons This allows one to reduce the problem 



to a known universality class of three dimensional models. 

Should such an approach (reducing the infrared dynamics to that of a known universality 
class) work in the case of the phase transition at Nf = N cr in QED 3 1 We believe that 
the answer is "no". The point is that there are too many relevant degrees of freedom at 
Nf = N CI in QED 3 . Besides a composite field <p, describing 2Nj Nambu-Goldstone (NG) 
bosons (corresponding to the spontaneous breakdown U(2Nf) — > U(Nf) x U(Nf)) and their 
2Nj — 1 flavor partners (for details, see below), there are massless fermions and massless 
photon at Nf = N cr . Therefore, unlike the temperature phase transition in QCD4, the 
effective action includes even more degrees of freedom than the initial lagrangian. This 
in turn implies that the effective action for the field (p, that can be obtained when light 
fermion and photon fields are integrated out, is necessary non-local and non-renormalizable. 
In other words, the universality class of the phase transition in QED 3 is different from those 
described by renormalizable ( or even local ) lagrangians for the field ip. 

Still, it would be interesting to find the effective action for <p in QED 3 , describing the 
phase digram of the theory. In this paper, we make a step in realizing this program and 
calculate the effective potential in leading order in 1/Nf expansion. Our main goal is to figure 
out how the structure of the potential reflects some peculiar properties of the phase transition 
at Nf = N CT in QED 3 pointed out in Refs. |^JT7|]. It gives an example of the conformal phase 
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transition |7j], the conception which generalizes the Berezinsky-Kosterlitz-Thouless (BKT) 
phase transition [18j ( taking place in two dimensions ) to higher dimensions. 

We stress that our aim is much more modest than to prove the existence of a finite 
N CT in QED 3 . We will just derive the effective potential, accepting, following Ref. ||, that 
the leading order in 1/Nf is a reliable approximation in the theory. Then we will see that 
the structure of this potential indeed corresponds to the conformal phase transition. In 
particular, it will be shown that unlike the usual, cx-model-like, continuous phase transition, 
there is an abrupt change of the spectrum of light excitations at the critical point in QED%. 
Also, the realization of the conformal symmetry is very different in the symmetric and 
non-symmetric phases. 



II. THE MODEL. GENERAL PROPERTIES. 



The lagrangian density for massless QED 3 is 

1 



€ = WW, (1) 

where = — ieA^, D = 7 M -D M , and four-component fermion fields carry the flavor index 
i = 1, 2, ... , Nf. The three 4x4 7- matrices can be taken to be 

^(?-° 3 )< (2) 

Recall that (cr 3 , iai, ia 2 ) and (— cr 3 , —ia\, — zcr 2 ) make two inequivalent representations of the 
Clifford algebra in 2 + 1 dimensions. 
There are two matrices 

that anticommute with 7°,7 1 and j 2 . Therefore for each four-component spinor, there is a 
global U(2) symmetry with generators 

I, \l\ 7 5 , and ^[ 7 3 ,7 5 ], (4) 

and the full symmetry is then U(2Nf). 

The lagrangian density with a mass term rmpip is invariant under parity transformations 
defined as 

P : ip(x ,x\x 2 ) -> t7 3 7V(z°, ~x\x 2 ). (5) 

The generation of a parity-invariant dynamical mass for fermions leads to the spontaneous 
breakdown of the U(2Nf) down to U(Nf) x U(Nf) with the generators 

£ (6) 
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a = 0, 1, ... , Nf — 1. The corresponding 2Nj NG bosons are 



(iV? scalars), and 



s « ~ ^ 7 5 V> (7) 



p« ~ ^ y7 3 ^ (8) 



(iV| pseudoscalars) . There are also their 2Nj massive flavor partners: 



(Nf scalars), and 



S a ~ ^ $ (9) 



P a ~^\[l\l^ (10) 



(NJ pseudoscalars). Notice that these 4Nj bosons can be described by the hermitian matrix 
field 

<f>) = ^Vj (11) 
which can be decomposed into a traceless part and its trace defined as: 

^0-^7°[7 5 ,7 3 k, X = M~7°[7 5 ,7 3 ^)~^ - (12) 

The field ip is assigned to the adjoint representation of SU(2Nf), and \ is a singlet. As to 
the vacuum group U(Nf) x U(Nf), the 2iV| NG bosons are assigned to the representation 
(Nf, Nf) © (iV?, A 7 /), and the 2Nj massive bosons are assigned to the representation (Nf x 
Nj, 1 ) © ( 1 , N f x Nf ) = (Nj - 1 , 1 ) © ( 1 , Nj - 1 ) © 2 ( 1 , 1 ) of its maximal semi-simple subgroup 
SU(Nf) x SU(Nf). One of the two singlets (1,1) corresponds to the pseudoscalar x = 
^[j 3 , / -f 5 }ip ~ P°, and another to the scalar a = ~ S 10 . The vacuum expectation value 
(0|er|0) is an order parameter describing the spontaneous breakdown U(2Nf) — > U(Nf) x 
U(N f ). 

Recall that besides the generation of a parity-invariant mass, corresponding to a c = 
(0|er|0) 7^ 0, there might be the generation of a {7(2Afy)-invariant mass corresponding to 
Xc = (0 1 x 1 0) 7^ 0, thus violating parity. Here we accept arguments of Ref. |Lj| in the 
support of a solution with a parity-invariant vacuum. 

In leading order in 1/Nf, in Landau gauge, the SD equation for the fermion propagator 
G(p) = (A(p)p — £(jo)) -1 is reduced to the equations @,[|: 



A(p) 



S(p)=m + 



7T 2 NfP 



kZ(k) k } kH(k) p 



dk— — — — — r — + dk 



k 2 + T?(k)p+^ J k 2 + Y?(k)k + 



(13) 
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where a = e 2 Nf, A is an ultraviolet cutoff, and, for generality, the bare fermion mass mo is 
introduced. 

Notice that the fermion mass function is neglected in the vacuum polarization in this 
approximation. The reliability of this approximation has been studied in Refs. 0,0. In 
particular, in Ref. ||15|| , to keep A(p)=l, a non-local gauge was used. It was shown then 
that including the mass function in the vacuum polarization does not lead to qualitative 
changes of the results of Ref. ||: the scaling law for the dynamical mass function near the 
critical point has the same form ( see Eq.(^) below), with a somewhat different value of 
N CT : N CT = 4.3 instead iV cr = 3.2. Moreover, this new value agrees with the result of the 
second paper of Ref. ||, incorporating the next-to-leading corrections in l/Nf to the SD 
equation in Landau gauge. 

That study suggests that neglecting a fermion mass function in the vacuum polarization 
is a reasonable approximation for a small (S « a) fermion mass, and this condition is 
fulfilled for Nf close to the critical point N CT . Therefore including the mass function in the 
vacuum polarization seems do not change the qualitative features of the phase transition at 
Nf = N cr . Since the main aim of this paper is to give a qualitative description of this phase 
transition, the use of this approximation seems appropriate. 

We will return to equation ( |T3"| ) in the next section. 



III. EFFECTIVE POTENTIAL FOR LOCAL COMPOSITE FIELDS IN QED 3 



We will consider the effective action (generating functional for proper vertices) of the 
local composite field <p in QED3. Actually we will calculate only the effective potential, 
i.e. the part of the action without derivatives of (p, though we will also describe properties 
of other terms in the effective action which are necessary for establishing the origin of the 
phase transition. 

In the derivation of the effective potential, we will closely follow Ref. ||20|| , where the 
effective potential for local composite fields in quenched QED^ was derived. 

The effective action is defined in the standard way. First, one introduces a generating 
functional for Green's functions of the field tp: 



Z(J) = exp(iW(J)) = J dr/exp \i / d 3 x{C{x) +tr [J(x)<p{x)]) 



(14) 



where the 77 integration is functional, J{x) is the source for <p(x), and C(x) is the lagrangian 
density (0) (the symbol 77 represents all the fundamental fields (fermion and photon ones) 
of the model). 

The effective action for the field <p is a Legendre transform of the functional W(J): 



T{p c ) = W{J) - / d 3 xtr [J{x) Vc {x)\ 



(15) 

where ip c (x) = (0|<£>(a;)|0). ^From Eqs.fll4]) and fllBD one finds that the following relations 
are satisfied: 

<Pc(x), (16) 



5J(x) 
5T 

5<p c (x) 



-J(x) 



(17) 
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The effective action T can be expanded in powers of derivatives of the field C : 



d 6 x 



-V{y c ) + -Z{ Vc )tr{d^ c d^ c ) + 



(18) 



where V(<p c ) is the efective potential. 

The calculation of the effective potential is reduced to finding the Legendre transform 
of the functional W(J) with the source J independent of coordinates x. The field ip c is a 
matrix in the flavor space describing AN 2 — 1 boson fields. 

Let us first consider the case when a constant source is introduced for the composite field 
a = tpip: C — > C + Jipip. Then the source term in Z(J) has the form of a bare mass term, 
with m = — J. Eq.(|i~7|) now becomes 



dV 
da. 



-mo. 



(19) 



Therefore 



V 



m (x)dx, 



(20) 



where m (cr c ) defines the dependence of m = — J on the condensate a c = (0|cr|0) = (O|^0|O), 
which is an order parameter of the spontaneous breakdown U (2Nf) —* U (Nf) x U(Nf). The 
function m (cr c ) can be defined from the SD equations fll3|). 
The effective potential V is calculated in the Appendix. It is: 



V 



A 2 N 2 f Hl 
32 



7 3 1, 
-- + — - - sin 29 



4is 2 



v 



4z/ 2 



cos 29 



+ 



A 2 JV>AEg(| + 



4tt 2 z/ 2 A 



sin 



(21) 



up 



where v = JX — 1/4, A = 8/n 2 Nf, A is an ultraviolet cutoff, A np = a/8 is an ultraviolet 
cutoff for nonperturbative dynamics, like A x = 4iiF n ~ lGev in QCD4 fl21~| , [22"! (for details, 
see the Appendix), Sq = H(p 2 )\ p= o is an infrared fermion mass parameter, 
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, A np e' 
v log — 



So 



<5 = 31og2 + --2, 



1 2v 
b\ = — arctan — , 



v 



and 



A 



2(1 + 



1T 2 Nf • 



T(l+iv) 



r(| + f)r(f + f) 



Notice that expression (pT|) for 1/ is valid at S << A np = a/8. 

The dependence of Sq on the condensate a c is defined from the equation: 



Or. 



(0|^|0) 



2NfAm, 







up 



7T- 



where the nonperturbative part of the condensate is 



(22) 



6 



np _ _J_ A 

4 np dp 



f a2 



p=A np 



^a 2 y ( d Fr 



iu 1 iv 3 p 2 



^0 



p=A n 



(23) 



F is a hypergeometric function (for details, see the Appendix). The parameter m is ex- 
pressed through E as 



AE 3 /2 J9/4 + p 2 



m 



A 



sin(z/log — + u(5 + 5i)). 



(24) 



Eqs. pi|) , (P^) , (P^) , and ( PH ) define as an implicit, and rather complicated, function of cr c . 
However, the phase diagram of the theory can be established by studying V as a function 
of the fermion mass Sq. It is convenient to rewrite Eq.(^) as 



A 2 A 2 £3 

V = f — - 

32 



77 + 



4z/ 2 



+ 1- 



4z/V 



\ 2 
cos(2fl + 2i/<Ji) - - sin(2fl + 2v5i 



^^(9/4 + 1^) . 2 
H A o o , sin (6 + udi). 



4vr 2 z/ 2 A 
Then the gap equation is 
dV _ A 2 A 2 S 2 



(25) 



up 



d£ 16 

4A(9/4 + z/ 



9 \ / 15 \ 

-7 + — J sin(# + v6i) + f 2f - —J cos(# + i/^ 



+ 



This equation yields the following solutions: 



[3sin(# + v8 x ) - 2z/cos(# + vSJ] \ sin(# + v5 x ) = 0. 



s[, re) =A np exp (-^ + ^o + ^i) 



, 71=1,2,..., 



(26) 



(27) 



si" ) =A np exp(-^ + 5 + 5 1 -^, 7i = l,2,... . (28) 

One can check that while all the solutions fl2"?D correspond to minima of V, solutions ( |2"8"D 
correspond to maxima of the potential. Actually, only the global minimum, corresponding 
to n = 1, defines the stable vacuum p2| . Therefore the dynamical mass is 



m dyn = S = A np exp ^~ + 5 + 6i \ . (29) 

As it has to be, it coincides with S of Ref . || , derived from the SD equation. 
Few comments are in order: 

1. Since expressions ([21] ) and (|25|) for V are valid only if S << A np , the solution (|29|) 
exists when < v « 1, i.e. when < n2 S Nf —\ « 1- Therefore, in this approxima- 
tion, there is a critical value Nj = N CT = ~ 3.24, and expression (|29| ) is valid in the 
near- critical, scaling, region. 
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2. The last term in expression (|2T|) for V is connected with a perturbative contribution, 
i.e. with ultraviolet dynamics at p » A np = a/8, where the dimensionless running 
coupling constant is weak (see Eq.(|32|) below). This term occurs because of the pres- 
ence of a source J = — m 7^ outside the extrema of V. On the other hand, since this 
term is proportional to sin 2 (6' + i/(5i), the dynamical mass fl29|), defined at the minimum, 
is independent of the perturbative contribution. Notice also that while this contribu- 
tion is proportional to Nf, the nonperturbative contribution in V is proportional to 
Nj. 

The perturbative term is divergent. It is connected with the point that we calculate 
the effective action for local composite operators, which is a generating functional for 
proper vertices (Green's functions) of these operators. It is known that the standard 
renormalizations do not remove divergences from such Green's functions (for example, 
see Sec. 12. 15 in the book |25|). These divergences are connected with perturbative 



short-range fluctuations, having nothing with infrared non-perturbative dynamics of 
bound states. The receipt of dealing with them in studying bound states (used, for 
example, in the case of QCD sum rules, dealing with Green's functions of colorless 
composite operators) is known: one has just to subtract them. Then we are led to the 
non-perturbative effective potential: 



A 2 N 2 f Y? n r 7 3 I . _ / 3 



32 



3 '4* , sin2 *-l 1 + 4^» COS2 * 



(30) 



It contains only contributions from the nonperturbative dynamics and is independent 
of the cutoff A. Notice that in the symmetric phase, with Nf > N cr , the parameters 9 
and v in Eq. (|30|) are imaginary. 

Henceforth the cutoff A is put equal infinity and we will consider the region of the 
nonperturbative dynamics with momenta p < A np . 

3. We have derived the effective potential for the composite operator a = ipip. In princi- 
ple, one can derive the effective potential for other 4Nj — 2 boson fields, described by 
the composite field <p, by introducing additional ANj — 2 sources independent of x The 
calculations of this potential at Nf > 1 are involved, and the expression for V should 
be very cumbersome. We have not succeeded in getting it. Fortunately, the effective 
action for a c , with all other fields taken equal zero, is sufficient for our purposes. The 
reason is the following. On the one hand, since a c is the only order parameter of the 
spontaneous breakdown U{2Nj) — > U(Nf) x U(Nf), the effective action for a c defines 
the mass of o particle both in symmetric and nonsymmetric phases. On the other 
hand, as was already pointed out in the previous section, in the symmetric phase, all 
the 4Nj — 1 bosons connected with the composite field (p are assigned to the same 
(adjoint) irreducible representation of the SU(2Nf). This point will be enough for us 
for proving that the phase transition in QED 3 is indeed conformal. 

IV. CONFORMAL PHASE TRANSITION IN QED 3 

The characteristic feature of the phase transition at Nf = N CT in QED 3 is that the 
scaling function 
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/(iV / ) = exp(-^ + 5o + 5i), u=U^L-l (31) 



(N CT — ^| in this approximation), has an essential singularity at Ay = iV cr . In particular, 
considering the scaling function f(z) as an analytic function of the complex variable z = Nf, 
one finds that while lim 2 ^ Zcr f(z) = as z goes to z CT from the side of nonsymmetric phase 
(Nf < N cr ), lim 2 ^ 2cr f(z) 7^ as z — > z cr from the side of the symmetric phase with 
Nf > N CI . This feature implies that the continuous phase transition in QED3 satisfies the 
criteria of the conformal phase transition (CPT) introduced in Ref. 0- ( Notice that since 
Nf appears analytically in the path integral of the theory, one can give a nonperturbative 
meaning to the theory with noninteger Nf ). 

As shown in Ref. M, the CPT is characterized by the following two general features: a) 
unlike the usual (cr-model-like) continuous phase transition, there is an abrupt change of the 
spectrum of light excitations at the critical point, though the CPT is a continuous phase 
transition; b) the realization of the conformal symmetry is very different in the symmetric 
and nonsymmetric phases: in particular, the equation for the divergence of the dilatation 
current, d^D^, has essentially different structures in those two phases. 

In this section, we will discuss how these two features are reflected in the effective action 
in QED 3 . 

QED3 is a super-renormalizable theory where the coupling constant is dimensional. 
However, as was pointed already in Ref. 0, in the infrared region, with p « A np = 
a/8, the conformal symmetry is a good symmetry in the symmetric phase with Nf > N CI 
and rridyn = (at least in leading order in l/Nf). The point is that in that region, the 
dimensional coupling constant a drops out from the SD equations for Green's functions: 
this, in particular, can be seen on the example of the SD equation for the mass function 
E(p) (see Eq. ( |A4j) in the Appendix). This point is also connected with the fact that the 
dimensionless running coupling constant 

ex 

a(p) = — ; (32) 

yFJ 8p(i + n(p)) v ; 

(n(p) is a polarization operator) has an infrared stable fixed point a = 1 in the symmetric 
phase @. 

Of course, in the nonsymmetric phase, because of the generation of the fermion dynamical 
mass, the conformal symmetry is broken. 

As it was already pointed out in the previous section, the parameter A np is an effective 
ultraviolet cutoff of the nonperturbative dynamics. Since, in leading order in 1/Nf, the 
polarization operator Tl(p) = a/8p = A np /p in the symmetric phase 0, the running coupling 
(|32|) is a = 1 at all values of p as A np — > 00. Therefore the dynamics is conformal invariant 
in that phase in this limit. 

The physical meaning of this limit is the following. The near-critical dynamics around 
Nf = N cr implies that, for such values of Nf, masses of light excitations have to be much less 
than the parameter A np . We are interested in the low energy dynamics of these excitations, 
when their energies are much less than A np . The question is whether there are such light 
excitations described by the composite field a (and, in general, by the field if). If they are, 
then the field Eq, describing their fluctuations around the minimum of the effective potential 
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V np ( see Eq. (|30|)), has also to be much less than A np . This in turn implies that the potential 
V np , describing the nonperturbative dynamics, has to be much less than A np . 

This leads us to studying the auxiliary "continuum" limit A np goes to infinity with E 
and being finite (the critical dynamics). As usual for critical dynamics ||22|| , such a limit 
can be interpreted as a renormalization of the parameter Nf. In other words, one has to 
find such a dependence of Nf on A np that as A np goes to infinity, the effective potential V np 
would be finite for finite So- 

It is not difficult to show from Eqs. fl29|) and (|30|) that in the nonsymmetric phase the 
non-perturbative effective potential is in this limit: 

^-^(4 + i io «| + s io « 2 D- < 33 > 



As it has to be, Nf, defined from Eq. (|29|), is fixed as A np — > oo: Nf = N cr . 

But what is the form of the potential at Nf = N CT from the side of the symmetric phase 
in this limit? Q It is easy to show that in the symmetric phase, at finite S , the potential V np 
(j30]) diverges as A np goes to infinity (on technical side, it is connected with the point that 
sin26 l and cos 26 in the potential become hyperbolic functions for Nf > N cr ). In particular, 
for Nf close to N cr , the potential is: 



A 2 N^ 
32 



10 A np e 5 ° , 3 1 2 A np e^° 

T" 21og ^T + 2 log 



oo (34) 



as A np — > oo. 

On the other hand, the kinetic term and terms with a larger number of derivatives (the 
structure of which is the same in the symmetric and nonsymmetric phases) have to be finite 
in this limit. Indeed, the most severe ultraviolet divergences always occur in an effective 
potential. Since the potential V np (|34T ) diverges only logarithmically, all other terms in the 
effective action are finite. 

This situation implies that, in the symmetric phase, there are no light particles described 
by the potential V np . If they exist, they are heavy (with M ~ A np ) and therefore decouple 
from infrared dynamics. 

In the present case, a boson is such a particle. However, since, as was already pointed 
out in Sec. 2, in the symmetric phase, all the 4Nj — 1 particles described by the field if are in 
the same (adjoint) representation of the SU(2Nf), all of them decouple there. Therefore, in 
the limit A np — > oo, the symmetric phase in massless QED 3 is a Coulomb-like, conformal- 
invariant phase, describing interactions between massless fermions and photons. 

On the other hand, the non-perturbative effective potential, and the whole effective 
action, for this composite field is finite as A np — > oo in the nonsymmetric phase (see Eq. (|33"D ) . 



1 Since Nf is dimensionless, it is actually a function of the ratio A np / [i, where /i is a renormalization 
group parameter of the dimension of mass. It can be introduced in different ways. For example, one 
condition might be that the potential V np is equal to (fi) s at So = /i. In the end of the calculations, 
if this limit exists, [i can be expressed through such physical parameters as the mass of a boson (in 
both symmetric and asymmetric phases) or the dynamical fermion mass in the asymmetric phase 
(see equation 
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This reflects the point that in that phase, besides photons and fermions, there are other 
light (with a mass M « A np ) particles: 2N"j massless composite NG bosons and 2Nj — 1 
massive composite bosons with a mass M ~ m^ yn = S << A np (in the scaling, near-critical, 
region). Moreover, the conformal symmetry is explicitly broken in the continuum limit in 
the nonsymmetric phase. Indeed, as follows from Eq.(K3), 



(O|0"D M |O) = (0|^|0) = 4K P (S C 



A 2 N 



Eq=So 



6 



(35) 



where 8% is the energy-momentum tensor. Therefore there is a conformal anomaly in this 
phase, and the vacuum expectation value (0|6^|0) plays here the same role as the gluon 
condensate in QCD A . 

Thus, in leading order in 1/Nj, the phase transition in QED% possesses all the char- 
acteristic features of the conformal phase transition. In particular, in agreement with the 
general conclusion of Ref. |7]], there is no Ginzburg-Landau like effective action describing 
the near-critical dynamics in both symmetric and broken phases. 

It is instructive to compare the dynamics in QED% with the dynamics in QED 3 with a 
Chern-Simons (CS) term. As is known [p3|j24l1 , the CS term enforces the phase transition to 
be first order. Therefore an abrupt change of the spectrum of light excitations at the critical 
point is natural in that case ||24|| . In QED% without the CS term, there is still an abrupt 
change of the spectrum, though the phase transition is continuous. 

It is amazing how closely the dynamics in QED 3 resembles the dynamics in quenched 
QED^ [p^ , p5[] . The scaling law for md yn in QED4 has the form (^) with A np replaced by 
cutoff A and N CT /Nf replaced by /a^\ where the dimensionless critical coupling a£9 ~ 1. 
The effective potential Q3"0| ) also closely resembles the effective potential calculated in Refs. 
20|,|7|j. This implies that long-range dynamics, provided by strong Coulomb-like (~ 1/r) 



forces, are essentially the same in these two models. 

Recently, the existence of the conformal phase transition (CPT) in quenched QED^ 
has been confirmed by calculating the spectrum of composites directly from Schwinger- 
Dyson equations for Green's functions of local composite operators ||26|| . It would be worth 
considering such an approach in QED 3 . 

The CPT is a conception extending the BKT phase transition [|T^] (taking place in 
two dimensions) to higher dimensions. To see this in QED 3 , it is convenient to consider the 
continuum limit A np — > 00 with md yn = ^0 being fixed as a renormalization of the parameter 



Then, the scaling law ([29|) implies the following (3 function at Nf < N a 



my 1 : 



Slog A 



up 




3/2 



(36) 



Clearly, this function has an ultraviolet stable fixed point at Nf = N CT . On the other 
hand, since at Nf > N CT , where m dyn = 0, the infrared dynamics with p << A np = a/8 is 
essentially independent of A np , the /3 function is identically zero. Therefore there is a line 
of fixed points at Nf > N CT . 

Such a renormalization group (RG) structure coincides with that corresponding to the 
BKT phase transition. 
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Of course, real QED3 has a fixed integer Nf and a finite A np = a/8. Still, the /3-function 
(j36|) is useful in the description of its dynamics. Let us consider a RG invariant quantity X. 
Then: 

dX ~ 9X W)#r = a (37, 



dlogA np <91ogA np PV / ^iV; 

This relation demonstrates that the (3 function ^(Nj 1 ) defines how X depends on the 
coupling constant a = 8A np . 

Thus the CPT yields extension of the BKT phase transition to higher dimensions: the 
crucial property of those phase transitions is the presence of an essential singularity in the 
mass (energy) gap at the critical point. However, there is an essential difference between the 
realization of the CPT in two and higher dimensions. While there cannot be spontaneous 
breakdown of any continuous symmetry in two dimensions [^7]], there is no such a restriction 
in higher dimensions. In particular, there is a genuine spontaneous flavor symmetry breaking 
in QED 3 . 



V. CONCLUSION 

In this paper we have analyzed the structure of the effective action in QED3 in leading 
order in 1/Nf and showed that it reflects the existence of the conformal phase transition at 
Nf = N CT . How much of this picture survives beyond the 1/Nf expansion is still an open 
issue, though results of the lattice computer simulations || are encouraging. 

The existence of such a phase transition, which is an extension of the BKT phase transi- 
tion to higher dimensions (2 + 1, in this case), is interesting in itself. Also, the CPT in QED% 
may be connected with nonperturbative dynamics in condensed matter, in particular, with 
dynamics of a non- Fermi liquid P8"|j2"9"|| . 

The most interesting issue is the connection of the phase transition in QED 3 with a 
possibility of the existence of an analogous phase transition in a (3 + l)-dimensional SU(N C ) 
gauge theory In QED S the critical number Nf = N cr separates two phases with 

very different dynamics. At Nf > N C1: there is a Coulomb phase describing interactions of 
massless photons and fermions; at Nf < N cr the rich dynamics with both spontaneous flavor 
symmetry breaking and confinement is realized. These two nonperturbative phenomena are 
intimately connected in this model: at Nf < N CI massive fermions decouple from infrared 
dynamics, thus leading to a potential growing with a distance, and therefore to confinement. 
As was already mentioned above, strong Coulomb-like forces play important role in providing 
the CPT in QED 3 . 

Is a similar picture realized in a (3 + l)-dimensional SU(N C ) gauge theory? What is the 
interplay between the dynamics provided by strong Coulomb-like forces and those connected 
with topologically nontrivial fluctuations, like instantons, monopoles, etc.? QED3, with its 
rich dynamics, is a fruitful laboratory for studying those complicated issues. 
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APPENDIX A: DERIVATION OF THE EFFECTIVE POTENTIAL 



In this Appendix we derive expression ( |2l| ) for the effective potential V. 

First, we need to recall some properties of the solution of the SD equation (|I3|). This 
equation will be used to extract the dependance of the parameter m on the condensate a c . 

Differentiating the second of Eqs. ([13|) with respect to p, one finds that satisfies the 
differential equation 



d 

dp 



V 



(p + a/8) 2 dY, 



a 



7T 2 N f p 2 + E 2 (p) 



2p + a/8 dp 

with infrared (IRBC) and ultraviolet (UVBC) boundary conditions: 



(Al) 



^ dp 



= (IRBC) 

p=0 

p(p + a/8) dT, 



2p + a/8 dp 



+ S(p) 



p=A 



Acffi 
2 ~dp 



+ S(p) 



p=A 



m (UVBC). 



(A2) 
(A3) 



Following Ref. 
follows: 



], one can approximate the equation 
8 p 2 £ 



1|) by writing it in two regions as 



d_ 

dp 



2 dT, 
^ dp 



n 2 N f p 2 + Y% 



2 ' 



p < A 



a 



up 



(A4) 



d 
dp 



3 dE 
^ dp 



2a 



p > A 



np; 



(A5) 



where, in Eq. (|A4[) , E 2 (p) in the denominator is replaced by Eq = £(p)| p= o; in Eq. (|A5|) , at 
p > A np , S 2 (p) in the denominator is neglected. This is justified if £ 2 Q») << A np : as we 
will see, it is indeed correct in the near-critical region. 

The parameter A np = a/8 is an effective ultraviolet cutoff for nonperturbative dynamics 
in the model: at tjiq = 0, the dynamical mass function rapidly decreases at p > A np (see 
below). This parameter plays here the same role as the parameter A x = 47rF n ~ lGev in 

qcd 4 

Taking into account IRBC (|A"3|), the solution to Eq. (|A"4T) is expressed through a hyper- 
geometric function: 



(I iv I iv 3 p 2S 



(A6) 
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where v 



■K 2 Nt 



— j. At p » S , its asymptotics is j3~0 



where 



and 



£,3/2 

£(p) = A — ^— sin [ulog — + v5 
vjp V E 



T(l + iu) 



A = ^ 
2 



r(l + |)r(| + f; 



(A7) 



(A8) 



) = - arg 
v 



r(| + f)r(f + f) 



(A9) 



(notice that <5 = 3 log 2 + | — 2 as i> — > 0). 

Solutions of Eq. (|A~5|) can in principle be expressed through Bessel functions, but, for our 
purposes, it is sufficient to display them as a series in x = p/at: 



S = x s [ C + — 

x 



Then we find from Eq. 



\ TT z J\f p J p Z 



(A10) 



m = C 1 + 



1 a 
7i 2 N f A / 



(All) 



Notice that Eqs. ( |A10|) and ( |A1 1|) imply that at mo = the dynamical mass function rapidly 
(as 4f) decreases with increasing p in the region with p > A np . Therefore, as was already 
stated above, A np ~ a is indeed an ultraviolet cutoff for nonperturbative dynamics. 

Matching now solutions of equations (|A4p and ( |A5| ) at p = A np (i.e. equating the values 
of the functions and their derivatives), one gets: 



£(p) 



a 



P=A np V 7T iy f A np y 



A2~ ; 



(A12) 



dp P=A np 

Eqs. (|ATI|) , (|AT2|) , and ( [All ) imply 



2Co a Q „ a 2 

— ZL-2" 



7T 2 N f A 2 A 3 ' 

J np np 



(A13) 



1 



m 



1 + 



r 2 iV f 



2 dp, 



P=A„ 



(A14) 



(here we neglect a term of order ■?). Then, using this equation and Eq. 
the equation 



we arrive at 
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m 



AE 3 J 2 19 



2ujA 



where 



Mip 



.4 



j + v 2 sin (v log + u(8 + 8i)j , 



(A15) 



An 



1 + 



1 2u 
di = — arctan — . 
v 3 



At m = we find solutions for the dynamical mass (compare with Eq. 

£ n) = A np exp {-~ + 5o + Si), n = 1,2, .... 
Eq.® implies that V function of E is: 

d(7c. 



(A16) 



V(E ) = - J m (x)-^-dx. 



(A17) 



The function mo (So) is given by Eq. ([A15 ). Let us determine the function %s The conden- 
sate a r is: 



a c = (0$(x)ip{x)\0) = -ti(0\Til>(x)$(y)\0) 

1 



-itrG(x, x) 



(20 



p-E(p) 



P 2 + £ 2 (p) 



(A18) 



The nonperturbative part of the condensate corresponds to the integration region with 



a 



up 



2N f 



71- 



dp 



p 2 E(p) _ Nj A2 dZ 



p 2 + Ti 2 (p) 4 np dp 



p=A np 



(A19) 



Here the last equality follows directly from equation 

Because, as follows from Eqs.( |AlO| ) and (All), £(p — > oo) = mo, the perturbative part 
of the condensate is <r Pt = —^fm A. This and equation ( |A19|) lead to expression (p2|) for 
the condensate a c = cr^ p + <r Pt . 

Eqs.(p2|), (|A15|) and ( |A17| ) lead to expression ( f2l|) for the effective potential. Notice that 
the expression for the potential in the symmetric phase, with Nf > N CT , is an analytic 
continuation of the potential in the broken phase, with Nf < N CI . 
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